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LONG-TERM  GOALS 

The  long-term  goal  is  to  improve  our  understanding  of  heat  and  moisture  exchange  in  the  tropics 
through  direct  estimates  of  the  fluxes  and  their  related  mean  variables.  The  flux  of  heat  across  the 
coupled  boundary  layers  is  primarily  accomplished  by  small-scale  processes  that  are  parameterized  in 
numerical  models.  The  ultimate  goal  is  to  improve  the  Navy’s  predictive  capabilities  in  the  tropics 
through  an  improved  understanding  of  the  processes  driving  the  Madden- Julian  Oscillation  (MJO). 

OBJECTIVES 

An  objective  of  the  LASP  program  is  to  improve  our  understanding  of  heat  and  moisture  exchange  in 
the  tropics  through  direct  estimates  of  the  fluxes  and  their  related  mean  variables.  The  flux  of  heat 
across  the  coupled  boundary  layers  is  primarily  accomplished  by  small-scale  processes  that  are 
parameterized  in  numerical  models.  The  primary  objectives  of  this  research  are  to: 

•  Deploy  a  suite  of  unique  observational  assets  in  and  over  the  Indian  Ocean  to  investigate  the 
MJO. 

•  Use  these  measurements  to  improve  the  surface  flux  parameterization  for  latent  and  sensible 
heat  used  in  Navy  models  and  LASP-related  observational  process  studies. 

•  Quantify  the  role  of  air-sea  exchange  in  the  formation  and  propagation  of  the  active  phase  of 
the  MJO. 

We  are  collaborating  with  researchers  from  NCAR,  NOAA/ETL,  Oregon  State  University,  and  other 
institution  to  investigate  the  relationship  between  boundary  layer  structure  and  surface  forcing  during 
an  MJO  event.  This  is  being  accomplished  through  measurements  collected  from  a  research  vessel 
that  conducted  surveys  during  three  20-40  cruises  to  investigate  air-sea  interaction  during  periods 
when  conditions  are  favorable  for  MJO  formation  (Madden  and  Julian,  1994).  The  measurement  will 
include  surface  meteorological  and  atmospheric  vertical  structure  and  collaboration  with  numerical 
modelers  and  other  observational  components  of  the  program.  The  principle  hypothesis  of  this 
research  is  that  improved  observations  and  parameterizations  of  latent  and  sensible  heat  fluxes, 
which  is  a  primary  source  of  energy  for  these  convective  systems,  will  improve  our  ability  to 
simulate  and  predict  the  MJO. 
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APPROACH 


Flux  Measurements:  The  PI  (Edson)  deployed  a  Direct  Covariance  Flux  System  (DCFS)  aboard  the 
R/V  Revelle  alongside  a  suite  of  instruments  to  measure  the  short  and  longwave  radiative  fluxes,  wind 
speed  and  direction,  temperature,  pressure,  humidity,  and  precipitation.  The  DCFS  (Edson  et  ah, 
1998)  has  been  used  in  a  number  of  field  programs  and  would  provide  estimates  of  the  momentum, 
sensible  heat  and  latent  heat  fluxes  during  ship-based  surveys.  The  PI  deployed  a  newly  developed  LI- 
COR  LI-7200  infrared  gas  analyzer  that  is  expected  to  improve  the  latent  heat  flux  estimates.  This 
new  unit  ran  alongside  the  LI-7500  that  has  been  successfully  deployed  in  previous  investigations. 

The  DCFS  and  infrared  hygrometers  allow  the  PI  and  his  colleagues  to  investigate  the  exchange  of 
sensible  and  latent  heat  between  the  atmosphere  and  ocean  using  the  direct  covariance  method.  This 
method  correlates  fluctuations  in  the  vertical  velocity,  W ,  with  fluctuations  in  the  sensible  heat, 
pac  T ,  and  latent  heat,  paLvq' ,  per  unit  volume: 


Qsen  =  PaCpW'T' 

(1) 

Qiat  =PaLy<l' 

(2) 

where  pa  is  the  density  of  air,  c  is  the  specific  heat  of  air,  Lv  is  the  latent  heat  of  vaporization,  V  and 
q'  are  temperature  and  specific  humidity  fluctuations,  respectively;  and  the  overbar  denotes  a  time 
average  ranging  between  10-30  minutes  for  turbulent  fluxes.  The  sensors  are  capable  of  accurately 
measuring  fluctuation  at  approximately  2  Hz  to  capture  the  total  flux  near  the  air-sea  interface. 

Unfortunately,  this  direct  method  is  generally  difficult  to  implement  over  the  ocean  due  to  platform 
motion,  flow  distortion  and  sensor  limitations.  Instead,  oceanographers  and  meteorologists  often  rely 
on  bulk  formula  such  as  the  COARE  3.0  algorithm  (Fairall  et  ah,  1996;  2003)  that  relates  the  fluxes  to 
more  easily  measured  averaged  wind  speed,  temperature  and  humidity.  These  averaged  variables  are 
related  to  the  flux  through  transfer  coefficients.  This  same  approach  is  commonly  used  to  parameterize 
the  surface  fluxes  in  forecast  models  from  variables  resolved  by  the  model.  For  example,  based  on  the 
dimensional  arguments,  the  exchange  of  sensible  and  latent  heat  at  the  ocean  surface  is  expected  to  go 
as  the  wind  speed  time  the  air-sea  temperature  and  humidity  differences,  respectively: 

Qsen  =  ~ PaCpCHUrbS 9  (3) 

Qiat  =  -paLvCEU,AQ  (4) 


where  CH  and  CE  are  the  transfer  coefficients  for  heat  and  mass  known  as  the  Stanton  and  Dalton 
numbers,  respectively;  U r  is  the  wind  speed  relative  to  water  (i.e.,  the  wind  speed-current  difference); 
and  A0  and  A Q  are  the  mean  air-sea  potential  temperature  and  specific  humidity.  The  uncertainty  in 
the  transfer  coefficients  for  heat  and  mass  remains  one  of  the  main  obstacles  to  accurate  numerical 
forecasts.  Improvement  of  these  transfer  coefficients  is  a  primary  objective  of  this  research. 
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WORK  COMPLETED 


Preliminary:  In  preparation  for  the  field  work  in  the  fall  of  201 1,  the  PI  combined  heat  flux  estimates 
from  a  number  of  recent  field  programs  to  look  at  the  behavior  of  the  transfer  coefficients  from  prior 
experiments.  These  field  programs  including  the  ONR  sponsored  CBLAST  program  (Edson  et  al. 
2007)  and  the  NSF  sponsored  CLIMODE  (Marshall  et  al.  2009)  and  GASEX  programs.  The  CBLAST- 
LOW  experiments  were  primarily  conducted  in  low  to  moderate  winds  while  the  CLIMODE  and 
GASEX  experiments  focused  on  air-sea  interactions  at  moderate  to  high  winds.  The  combined  data  set 
therefore  covers  a  wide  range  of  wind  and  stability  conditions.  For  example,  near  surface  winds  of  15 
m/s  were  commonly  encountered  over  the  North  Atlantic  during  CLIMODE  and  the  data  set  includes 
wind  events  with  speeds  over  25  m/s.  These  high  wind  events  drive  surface  stresses  that  routinely 
exceed  1.0  N/m  and  combined  latent  and  sensible  heat  fluxes  from  the  ocean  into  the  atmosphere  that 
exceed  1200  W/m2.  These  enormous  heat  fluxes  are  driven  by  high  winds  and  large  air-sea  temperature 
and  humidity  differences  encountered  over  the  Gulf  Stream  during  cold  air  outbreaks.  The  CBLAST- 
LOW  program  collected  3  months  of  data  from  an  Air-Sea  Interaction  Tower  (ASIT)  under  low  to 
moderate  wind  conditions.  To  data,  the  CBLAST  investigations  have  focused  on  the  role  of  swell  on 
momentum  exchange  under  low  wind  conditions  while  the  CLIMODE  investigations  have  focused  on 
momentum  exchange  at  high  winds.  This  investigation  focuses  on  heat  exchange  from  low  to  high 
winds  using  the  combined  data  sets.  Preliminary  results  from  this  investigation  were  reported  at  the 
AMS  17th  Conference  on  Air-Sea  Interaction  in  Annapolis,  MD. 

Main  Experiment:  The  PI  (Edson)  deployed  a  turbulent  and  radiative  flux  package  and  associated 
mean  meteorological  sensors  aboard  the  research  vessel  the  R/V  Revelle  during  the  DYNAMO  field 
program.  In  situ  meteorology  and  high-rate  flux  sensors  operated  continuously  while  in  the  sampling 
period  for  DYNAMO  Leg  3.  This  included  all  sensors  operating  during  Leg  2  with  the  addition  of  a 
closed-path  LI-7200  IRGA  to  the  flux  systems.  Sea  surface  temperatures  were  measured  by  the  group 
using  the  sea-snake  floating  thermistor  and  radiometric  estimates  of  skin  temperature  in  collaboration 
with  Chris  Zappa  (LDEO/Columbia).  NOAA/PSD  operated  a  suite  of  remote  sensing  instruments  for 
low  clouds  and  light  precipitation:  the  NOAA  W-band  cloud  radar,  a  microwave  radiometer,  and  a 
laser  ceilometer.  Aircraft  overflights  were  made  on  November  13,  22  and  26  with  all  systems 
operational  and  good  relative  winds  for  our  flux  measurement  systems. 

Overall,  these  packages  ran  continuously  in  international  waters  from  4  September  (start  of  Leg  1)  to 
3 1  December,  2011  (end  of  leg  4).  The  data  return  rate  was  nearly  100%.  Preliminary  results  from 
the  field  campaign  were  reported  at  the  AMS  18th  Conference  on  Air-Sea  Interaction  in  Boston,  MA. 

The  data  has  been  quality  controlled  and  calibrated  to  produce  means  and  bulk  fluxes  using  the 
COARE  3.5  algorithm  (Edson  et  al.,  JPO,  43,  1589-1610,  2013).  These  data  are  available  through  an 
anonymous  ftp  site  maintained  at  UConn:  ftp://dynamo.dms.uconn.edu/.  These  data  include  radiative 
and  heat  fluxes;  surface  stress;  wave  data;  surface  and  near  surface  sea  temperatures,  salinity  and 
currents;  and  other  key  variables  specifically  requested  by  DYNAMO/LASP  Pis. 

RESULTS 

Preliminary  Research:  The  observational  highlight  of  Leg  3  was  the  capture  of  an  almost  complete 
MJO  cycle  in  our  time  series  measurements.  After  a  brief  period  of  convective  activity  upon  arrival, 
the  measurements  were  characteristic  of  the  suppressed  phase  of  the  MJO.  Strong  solar  heating  of  the 
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ocean  overwhelmed  the  convective  cooling  by  the  atmosphere,  and  ocean  surface  temperature 
increased  by  approximately  1°C  between  1 1-17  November  (Yeardays  315-321).  As  shown  in  Figure 
3,  winds  remained  light  and  variable  during  this  period  and  little  precipitation  was  observed.  In 
Figure  3,  the  velocity  data  were  taken  from  the  sonic  anemometers  on  the  forward  mast.  Poor  relative 
wind  directions  were  removed  and  the  data  was  interpolated  through  these  periods.  Surface  currents 
were  collected  from  the  ship’s  ADCP  after  QC  by  the  OSU  mixing  group.  Air  temperature  was 
collected  by  two  aspirated  T/RH  sensors  on  the  bow.  Poor  relative  wind  directions  were  also  removed 
from  these  data  to  reduce  the  heat  island  effect  of  the  ship.  Sea  surface  temperature  was  measured  by 
the  sea-snake  and  corrected  for  cool-skin  effects.  Precipitation  (P)  was  provided  by  an  optical  rain 
gauge  after  calibration  with  6  other  rain  gauges  (see  below).  Specific  humidity  was  calculated  by 
combining  these  sensors  with  pressure  measurements.  Evaporation  (E)  was  computed  using  the  latent 
heat  flux  estimated  by  the  COARE  3.5  algorithm.  Accumulated  totals  of  P  and  E  are  shown.  The 
blue  shaded  area  is  therefore  their  difference. 
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Figure  3.  Time  series  of  mean  meteorological  and  ocean  surface  variable  during  Leg  3. 


Atmospheric  convection  began  to  increase  on  17  November  (Yearday  321)  with  precipitation  falling 
overnight.  Wind  speeds  and  convection  gradually  increased  during  the  period  between  18-23 
November  (Yeardays  322-327).  Sea  surface  temperatures  leveled  off  during  this  period  due  to  a 
gradual  reduction  in  solar  radiation  and  continued  latent  heat  loss  of  approximately  100  Wm'2.  The 
active  phase  of  the  MJO  was  in  full  swing  with  the  arrival  of  cyclone  aided  westerly  wind  burst  on  24 
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November  (Yearday  328).  Ten-minute  averaged  wind  speeds  in  excess  of  19  m/s  where  observed 
over  growing  seas  with  significant  wave  heights  of  approximately  3  m.  The  sea-surface  temperature 
dropped  significantly  during  this  period  with  combined  sensible  and  latent  heat  loss  to  the  atmosphere 
of  approximately  200-400  Wm'2.  Moderate  winds,  precipitating  convection  and  surface  cooling 
continued  through  the  end  of  November  (Y earday  334).  Overall,  precipitation  exceeded  evaporation 
during  this  period.  Although  westerly  surface  winds  around  5  ms'1  continued  into  December,  a 
drying  out  of  the  atmosphere  aloft  and  associated  reduction  in  convective  activity  was  observed  as  the 
active  phase  of  the  MJO-related  convection  moved  towards  the  maritime  continent. 

Surface  Fluxes:  A  summary  of  the  latent,  sensible,  radiative  and  net  heat  fluxes  are  shown  in  Figure 
4.  The  downward  radiative  fluxes  were  measured  by  the  purgeometers  (LW)  and  pyranometers  (SW) 
located  on  the  top  of  the  forward  mast.  The  upward  long-wave  radiation  was  obtained  using  our  SST 
measurements  and  corrected  for  IR  reflection  from  downwelling  long-wave.  A  commonly  used 
parameterization  of  sea-surface  albedo  was  used  to  estimate  the  reflected  solar.  The  optimized  set  of 
mean  meteorological  and  surface  ocean  measurements  (temperature  and  currents)  are  used  to  compute 
the  latent,  sensible  and  rain  fluxes  with  the  CO  ARE  3.5  algorithm. 


I 


500 

0 

-500 

■1000 


315  320  325  330  335  340 


xIO7 


Figure  4.  Time  series  of  radiative,  surface  heat  and  net  heat  fluxes  during  Leg  3.  The  lower  panel 
shows  the  integral  of  the  net  heat  flux  over  the  experiment.  A  small  correction  (~1  %  reduction) 
was  applied  to  the  short-wave  measurements  from  the  forward  mast  based  on  comparison  with  the 
PSD  and  OSU  sensors.  A  more  significant  correction  to  the  long-wave  measurements  was  applied 
to  the  mast  using  the  PSD  sensors.  The  latent,  sensible  and  rain  induced  fluxes  are  computed 

using  CO  ARE  3.5. 
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The  lower  plot  in  Figure  4  shows  the  integrated  value  of  the  net  heat  flux  into  the  ocean.  The  large 
amount  of  heat  supplied  to  the  ocean  during  the  suppressed  phase  of  the  MJO  is  clearly  seen  in  the 
measurements  between  1 1-20  November  (Yeardays  315  and  324).  It  is  important  to  note  that  not  all 
of  this  energy  goes  into  heating  of  the  mixed  layer  as  a  significant  fraction  of  the  solar  radiation 
penetrates  through  this  layer.  Oceanic  turbulence  also  removes  some  of  the  heat  through  the  base  of 
the  mixed  layer.  However,  this  still  represents  a  significant  amount  of  energy  stored  by  the  oceanic 
capacitor  to  potentially  initiate  (in  tandem  with  favorable  atmospheric  condition,  e.g.,  surface 
convergence  and  divergence  aloft)  and  drive  convection  during  the  active  phase  of  the  MJO.  This  plot 
also  shows  a  leveling  off  of  the  heating  around  21  November  (Yearday  325),  which  likely  corresponds 
to  the  start  of  the  active  phase.  The  westerly  wind  burst  is  clearly  associated  with  the  active  phase, 
and  drove  the  net  heat  loss  evident  between  24  November  and  1  December  (Yeardays  327  and  335) 
and  associated  cooling  of  the  ocean.  The  end  of  the  time  series  again  shows  the  heat  supplied  to  the 
ocean  (i.e.,  recharging  of  the  capacitor)  during  the  start  of  a  new  suppressed  phase.  The  rate  of 
energy  input  is  very  similar  to  the  prior  period  and  represents  the  clear  sky  solar  value  integrated  over  a 
day  minus  the  nominal  value  of  latent  heat  flux  of  approximately  100  Wm"2.  The  role  that  the  ocean 
and  air-sea  exchange  plays  in  driving  the  amplitude  and  phase  of  the  MJO  will  be  a  focus  of  our 
collaborative  investigation  with  the  Ocean  Mixing,  Sounding  and  Remote  Sensing  groups. 
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Figure  5.  Time  series  of  surface  stress  (top)  and  latent  heat  flux  (bottom)  derived  from  the 
direct  covariance  and  bulk  aerodynamic  methods.  The  solid  color  is  the  bulk  and  the  black  line 

is  the  direct  measurements. 

Surface  Flux  Parameterization:  The  flux  group  has  begun  to  process  the  turbulence  instrumentation 
on  the  forward  mast  to  compute  fluxes  using  the  direct  covariance  (eddy  correlation)  method  after 
correction  for  ship  motion.  Preliminary  results  from  our  analysis  are  encouraging.  For  example, 
Figure  5  shows  a  time  series  of  surface  stress  and  latent  heat  flux  estimated  from  the  direct  covariance 
method  and  COARE  3.5  bulk  algorithm  during  the  westerly  wind  burst.  In  this  figure,  the  solar 
background  color  represents  the  bulk  estimates,  while  the  solid  line  represents  the  DC  measurements. 
The  DC  measurements  are  limited  to  more  favorable  relative  wind  directions  to  result  the  impact  of 
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flow  distortion  on  our  measurement.  Nonetheless,  the  agreement  between  the  two  stress  estimates  is 
very  good. 

A  primary  goal  of  our  research  is  to  improve  the  surface  flux  parameterization  for  latent  and  sensible 
heat  used  in  these  models  and  observational  process  studies.  The  starting  point  of  this  investigation  is 
the  comparison  of  the  directly  measured  transfer  coefficients  of  sensible  heat  and  moisture  known  as 
the  Stanton  and  Dalton  numbers,  with  the  COARE  3.5  algorithm  (Edson  et  al.  2013).  In  CO  ARE  3.5 
the  coefficients  of  sensible  heat  and  moisture  are  assumed  to  be  the  same,  suggesting  similarity  in  the 
transfer  of  heat  and  mass.  The  investigation  will  test  the  assumption  of  similarity  and  the  consistency 
of  these  results  with  previous  estimates  of  the  Stanton  and  Dalton  numbers  from  the  CBLAST, 
CLIMODE  and  GASEX  programs. 

An  example  of  this  effort  is  shown  in  Figure  6  where  the  transfer  coefficients  computed  from  DC 
latent  heat  and  buoyancy  fluxes  during  DYNAMO  and  CLIMODE  are  compared.  The  four  estimates 
of  the  transfer  coefficients  are  in  very  good  agreement.  All  four  estimate  exhibit  the  same  overall 
characteristics  of  a  minimum  around  3-4  m/s  that  increases  to  a  maximum  in  the  transfer  around  12 
m/s  before  again  falling  off  at  higher  winds.  The  ultimate  goal  of  this  research  is  the  development  of 
the  COARE  4.0  algorithm  and  a  qualitatively  determined  parameterization  that  captures  these 
characteristics  is  drawn  on  this  figure.  A  rigorous  evaluation  of  this  function  is  being  conducted  using 
these  results  combined  with  those  from  CBLAST  and  GASEX. 
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Figure  6.  Transfer  coefficients  for  moisture  and  buoyancy  computed  from  DC  fluxes  and  their 
associated  means.  The  two  sets  of  transfer  coefficients  were  generated  from  data  collected  during 

DYNAMO/LASP  and  the  CLIMODE  programs. 
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The  improved  flux  estimates  are  being  used  to  investigate  boundary  layer  response  to  air-sea 
interaction.  For  example,  the  exchange  of  heat  across  the  air-sea  boundary  layers  plays  an  important 
role  in  the  thermodynamics  of  the  coupled  ocean-atmospheric  system  in  the  tropics.  Of  particular 
interest  is  the  relative  magnitude  of  the  net  shortwave  (solar)  radiative  flux  that  heat  the  upper  ocean 
versus  the  latent  heat  flux  that  cools  it  and  drives  atmospheric  convection.  All  the  terms  in  the  surface 
heat  budget  for  Leg  3  are  shown  in  Figure  7. 

The  net  heat  flux  shows  clear  phasing  with  the  active  and  inactive  phases  of  the  MJO  with  a  release  of 
energy  during  the  active  phase  and  oceanic  storage  of  energy  between  events.  For  example,  the  low 
panel  shows  the  integrated  value  of  the  net  heat  flux  into  the  ocean  in  J/m2.  One  can  think  of  the  ocean 
as  a  capacitor  for  the  MJO,  where  the  energy  is  being  accumulated  when  there  is  a  net  heat  flux  into 
the  ocean  (here  occurring  to  approximately  November  24)  after  which  it  is  released  to  the  atmosphere 
during  the  active  phase  of  the  MJO  under  high  winds  and  large  latent  heat  exchange. 


Figure  7.  SST  (upper  panel),  heat  and  radiative  flux  (2nd  panel),  net  heat  flux  (3rd  panel)  and 
integrated  energy  flux  (lower  panel)  from  data  collected  during  Leg  3. 

The  active  phase  of  the  MJO  can  also  drive  strong  wind  events  and  oceanic  response  that  include  rapid 
enhancement  of  the  equatorial  jet  and  deepening  of  the  mixed  layer.  Interestingly,  the  SST  (upper 
panel)  is  decreasing  prior  to  the  reversal  of  the  net  heat  exchange  out  of  the  ocean.  The  physical 
processes  responsible  for  this  decrease  are  being  explored  with  our  collaborators  at  OSU  and 
NOAA/PSD.  An  overview  paper  and  a  modeling/observation  study  both  describing  our  preliminary 
investigations  of  the  MJO  disturbance  captured  during  Leg  3  have  been  submitted  to  BAMS  (Mourn  et 
ah,  submitted  summer  2013;  Chen  et  ah,  submitted  summer  2013) 
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IMPACT/APPLICATIONS 


None  to  date 

TRANSITIONS 

None  to  date 

RELATED  PROJECTS 

The  ONR  portion  of  this  program  will  work  closely  with  investigators  from  the  NSF/NOAA  funded 

DYNAMO  program.  The  results  will  also  be  compared  with  findings  from  the  NASA/SPURS 

program. 
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